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Abstract 


We used the light curve code XRBinary to model the quiescent K2 light curves of three low-inclination cataclysmic 
variables (CVs): IRXS J06324-2536 (J0632--2536), RZ Leo, TW Vir and the pre-CV WD 1144+011. Optimized 
light curve models were obtained using a nonlinear fitting code NMfit and visualized by Phoebe 2.0. The disk 
model of J06324-2536 shows that one hotspot at the edge of the disk is enough to describe its light curve, while the 
other two dwarf nova (DN): RZ Leo and TW Vir require two hotspots. A typical pre-CV model with a weak 
irradiation effect for WD 1144+011 can explain its single-hump modulation and the newly observed spectrum 
confirms its previous classification. The synthetic analyses for the DN clearly indicate that phase zero of the 
double-hump modulations occurs around the secondary minimum and the primary hump is mainly caused by the 
hotspot at the edge of the disk. The quiescent disk has a flat temperature distribution with a power index of ~0.11. 
The disk model of RZ Leo implies a truncated disk, supporting its previously speculated classification as an 
intermediate polar (IP). Except for the IP model of RZ Leo, which lacks a component related to the inferred 
accretion curtain, the models of J06324-2536, TW Vir and WD 11444-011 are consistent with results from the Gaia 
mission. The derived masses and radii of the secondaries of the three DN are consistent with the semi-empirical 
relations for CV donor stars, while their effective temperatures are higher than the predictions. Irradiation of the 
donor stars is investigated to explain this discrepancy. 


Key words: binaries (including multiple): close — novae, cataclysmic variables — white dwarfs 


1. Introduction 


Dwarf novae (DN) are a subtype of primarily non-magnetic 
cataclysmic variables (CVs), in which a white dwarf primary 
accretes matters from a Roche-lobe filling late-type star via the 
inner (L1) Lagrange point (Warner 2003). An accretion disk can 
extend to the white dwarf surface via viscous processes between 
adjacent accretion annuli (e.g., friction and shear) due to the 
weak magnetic field of the white dwarf in DN (B « 106 G). 
Systems in which the magnetic field of the white dwarf is large 
enough (10° « B « 107 G) for the accretion disk to be disrupted 
inside the white dwarf magnetosphere and for material to begin 
to follow the magnetic field lines are called Intermediate Polars 
(IPs). In the following, we assume that the disk outside of the 
magnetosphere in an IP is equivalent in structure to the disk in a 
regular non-magnetic DN. The interaction between the ballistic 
stream leaving the L1 point and the accretion disk forms a region 
of energy release at the edge of the disk, which is called a 
hotspot. The combination of an accretion disk and a hotspot 
(hereafter just called the disk model) has been used as the typical 
accretion model for CVs (e.g., Smak 1970; Warner 2003) and 
has been quite successful in describing the asymmetrical eclipse 
profile of several high-inclination CVs during quiescence (e.g., 


Wood & Crawford 1986; Bruch 1996). This model can also be 
used to explain a symmetrical CV eclipse profiles during 
outburst (e.g., Kato et al. 2003; Bakowska & Olech 2015), 
which is caused by the significant flux increase of the disk (1.e., 
the luminous accretion disk almost overwhelms the relatively 
faint white dwarf and hotspot). 

Several high-inclination CVs in quiescence have been 
comprehensively studied using this disk model (e.g., Cook & 
Waner 1984; Bailey & Cropper 1991; McAllister et al. 2015). 
By detecting variations in the mid-eclipse times, substellar 
objects have been suggested to exist in several DN (e.g., 
V2051 Oph, Qian et al. 2015; and EM Cyg, Dai & Qian 20102). 
In addition, many synthetic light curve analysis methods (e.g., 
BINSYN program, Linnell et al. 2012; Eclipsing Light Curve 
Code, Orosz & Hauschildt 2000; and the cool-disk model, 
Khruzina 2011) have been developed to analyze CV eclipse 
light curves. The irregular eclipse light curves of quiescent CVs 
are composed of the occultations of multiple components, 
including the white dwarf, the accretion disk and hotspots (e.g., 
Smak 1994; Feline et al. 2004; Littlefair et al. 2014; McAllister 
et al. 2015). High-inclination CVs only show a single eclipse 
in one orbit (ie. the secondary eclipse of the red dwarf is 
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invisible; e.g., Krzeminski 1965; Bailey et al. 1988) because 
the red dwarf in a CV system is usually regarded to be a very 
faint component compared with the accretion disk and the 
white dwarf. This single-eclipse feature implies that CV 
eclipse light curves cannot reveal full information about the 
red dwarf companion. Due to the complexity of the quiescent 
CV eclipse light curves, high time resolution is necessary to 
decompose all components. However, over 70% of CVs with 
short orbital periods (<3hr) are fainter than 17 mag in 
quiescence based on the updated CV catalog (RKcat Edition 
7.24, which was first published in Ritter & Kolb 2003). 
Therefore, the majority of CVs have only low time-resolution 
photometry comprised of blended flux from many components, 
which makes modeling these systems difficult. For example, 
recent discussions concerning two new-found CVs with deep 
eclipses carried out by Kjurkchieva et al. (2015) and Kennedy 
et al. (2016) have clearly indicated that the model light curves 
cannot perfectly fit the observed eclipsing light curves during 
ingress and outside of eclipse. Assuming that a deep eclipse is a 
common feature of CVs with orbital inclination higher than 80°, 
the fraction of low-inclination CVs can be simply estimated to be 
around 90%. To determine a general model independent of 
inclination, it is necessary to consider low-inclination CVs. 

The unprecedented light curves from the Kepler K2 mission 
(Howell et al. 2014), with nearly continuous photometric 
coverage for 1-3 months at different pointings (Campaigns) 
along the ecliptic provide an excellent database to study 
quiescent CV light curves. K2 Campaign 0 (K2-CO) was an 
engineering test in the early stage of the K2 program and only 
covered ~35 days because the spacecraft was not in fine point 
during the beginning of the campaign, while K2 Campaign 
1 (K2-C1) covered a complete period of 80 days. In this paper, 
we focus on the phased light curves of four systems: 
1RXS J0632+2536 (hereafter J06324-2536) and TW Vir, 
which are both DNs; RZ Leo, which is an unusual IP that 
also has displayed DN outbursts, making it one of the few 
systems to be both a magnetic system and a DN; and WD 1144 
+011, which is a pre-CV (meaning that the secondary is 
probably not filling its Roche lobe). J06324-2536 was observed 
in K2-CO, while TW Vir, RZ Leo and WD 1144+011 were 
observed in K2-C1. The phased light curves are investigated in 
detail using the synthesis methods XRBinary and NMfit 
described in Section 3. Due to the lack of any eclipse feature, 
all four objects are likely low-inclination systems. The preset 
model parameters are discussed in Sections 4 and 5. The details 
of the white dwarf accretion structure during quiescence and 
the physical parameters of the stars in each system are further 
discussed and visualized in Section 6. 


2. Phased Light Curves 


Dai et al. (2016) used the PyKE suite of software tools 
developed by the Guest Observer Office (Still & Barclay 2012) to 
extract the K2 light curves of J0632--2536, RZLeo and 
WD 11444-011, and derive their orbital periods using traditional 
period finding techniques (e.g., Lomb-Scargle periodogram, 
Lomb 1976; Scargle 1982; and phase dispersion minimization, 
Stellingwerf 1978) and their corresponding phased light curves. 
Due to the unstable orbital modulations of TW Vir, Dai et al. 
(2017) applied a phase-correcting method to obtain its orbital 
period and the phased data from the quiescent data. 

Based on the continuous K2 data, we attempted to model the 
mean orbital light curves of each system by using an XRBinary 
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light curve synthesis code that was developed by E. L. 
Robinson.’ The default phase zero of a light curve generated 
by XRBinary is an inferior conjunction of the Roche-lobe 
filling secondary (ie., the accretion disk around the white 
dwarf is occulted by the red component), which is in accord 
with the phases specified in eclipsing CVs. Compared with the 
typical CV eclipse light curve with a narrow and deep white 
dwarf eclipse (e.g., Littlefair et al. 2014; McAllister et al. 
2015), all four light curves derived by Dai et al. (2016, 2017) 
only show wide and shallow modulations with an amplitude of 
several hundredths up to tenths of a magnitude. Note that the 
phases of the obtained light curves are arbitrary. The three DN 
(06324-2536, RZ Leo and TW Vir) light curves clearly show 
double-hump modulations, with a nearly constant phase 
difference of ~0.5 between the two minima dips. This means 
that the light minima at the lower (primary dip) and higher 
(secondary dip) flux levels are at phases 0.5 or zero, 
respectively. In principle, an irradiation effect is significant in 
CV systems and commonly results in a higher flux level at the 
phase 0.5 dip than at phase zero. However, phase zero of the 
double-hump modulation cannot be simply identified from this 
flux difference due to possible changes caused by hotspots on 
the disk. Therefore, for each DN with a double-hump 
modulation discussed in this paper, two phased light curves 
corresponding to phase zero at the primary and secondary dips, 
respectively, were analyzed to search for a final convergent 
solution. Since the secondary dip in TW Vir cannot be 
accurately measured due to large scatter in the data, the 
primary maximum was set to be at the reference phases 0.75 
and 0.25, guaranteeing that minima would occur around phases 
0.5 and zero. For the pre-CV WD 11444-011, phase zero was 
set to be the minimum of the light curve since the light curve 
shows a single-hump modulation. 

The four K2 light curves are expressed in Simple Aperture 
Photometry flux (i.e., electrons per second), while the light 
curves calculated by XRBinary are given in ergs per second. 
Hence, the K2 light curves were normalized before the 
calculations of XRBinary. Since the huge number of data 
points in the observed light curves (i.e., the number of data 
points in the K2 light curves observed in long cadence (LC; 
30 minutes sampling) and in short cadence (SC; 1 minutes 
sampling) are more than 10° and 10°, respectively) can take a 
long calculation time with XRBinary, the K2 phased light 
curves are moderately binned with a uniform phase resolution 
0.01. Dai et al. (2017) demonstrated that the orbital modula- 
tions of the three CVs are stable, and the large-amplitude 
dispersion of ~0.43 mag in the phased light curve of RZ Leo 
(Figure 5 of Dai et al. 2016) is only caused by a uniform drift of 
the system light. Thus, while the phased light curve of TW Vir 
has the highest stability (Dai et al. 2017), all four binned and 
normalized K2 light curves represent their orbital modulations. 


3. Synthesis Methods for the Light Curves 


Although the XRBinary program was initially designed to 
calculate the light curves of low-mass X-ray binary stars 
(LMXBs) (Gomez et al. 2015), it can also be used to model 
low-inclination CVs based on the following two reasons. First, 
XRBinary sets the primary star as a much smaller sphere than 
all other dimensions in the binary system and is unresolved by 
any of the grids used in calculating the light curves. Like the 


12 http: / [www.as.utexas.edu/ -elr/Robinson/XRbinary.pdf 
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neutron star in LMXBs, the white dwarf in both CVs and pre- 
CVs is very small compared with the companion red dwarf. 
Second, XRBinary simply assumes that the primary star only 
emits black body radiation. This is a reasonable assumption for 
low-inclination CVs and pre-CVs because the flux contribution 
from the visible white dwarf is approximately constant. 
XRBinary may not be appropriate to accurately model 
complicated eclipse light curves but it should be able to 
reproduce low-inclination CV light curves. In fact, XRBinary is 
a powerful tool to analyze the accretion disk around the central 
compact star by constructing a complex accretion disk model 
consisting of a disk, a disk rim, a disk torus and an inner disk. 
Ratti et al. (2013) have successfully applied XRBinary to 
reproduce an ellipsoidal light curve of the CV CXOGBS 
J174444.7-260330 in a low state. 

In this paper, a standard CV model (i.e., a semi-detached 
close binary system with a Roche-lobe filling red dwarf and, in 
the case of the 3 DN, an accretion disk around the white dwarf) 
with a set of geometric and physical parameters is applied to 
model the phased light curves of four low-inclination systems. 
Based on a complete set of preset parameters, XRBinary 
calculates a theoretical light curve and a relative measurement 
of goodness of fit, X? (i.e., the variance between the calculated 
and observed light curves). Taking advantage of the improved 
Nelder-Mead method (i.e., Multi-directional Search method 
(MDS), Nelder & Mead 1965; Torczon 1989; Geng 2008), 
which was successfully applied to fit the O-Cs of AM Her (Dai 
et al. 2013) and UZ For (Dai et al. 2010b), we developed a new 
program called NMfit to carry out a search for the best binary 
model within a given parameter space. After obtaining the best 
binary model, NMfit sets a series of small deviations around 
each of the parameters to test the variations in x^ caused by 
these deviations. Note that when an adjustable parameter is 
intentionally set to deviate from its optimal value, then the 
other parameters are fixed. Moreover, this optimization method 
ignores degeneracies between parameters, which is a big issue. 
Since the best binary model means the minimal x. the 
deviations give rise to an increase of x^. When X? increases to 
5% larger than the minimal value, this tested parameter 
deviation is arbitrarily regarded as the uncertainty. All of the 
adjustable parameters are tested one-by-one to estimate their 
errors. The error estimates based on x? are only lower limits to 
the true uncertainties in the derived parameters and the true 
uncertainties could be quite large because the physics of 
compact binaries is much more complicated and XRBinary 
arbitrarily makes many tacit assumptions (e.g., the disk is 
modeled as an object with sharp, well-defined edges and 
surfaces, and the disk surface elements emit and absorb like 
black bodies). The uniform design proposed by Fang (1980) is 
used to decide the initial parameter set before the search 
iterates. A good initial parameter set is crucial for the success of 
the iterations. In addition, since Phoebe? is able to calculate 
and plot an accretion disk around a compact star (Prša 
et al. 2016), it was used to visualize the system configuration 
based on the best-fitting parameters derived by NMfit and 
XRBinary. 


4. Model Parameters Preset in XRBinary 


The 23 model parameters that were used to generate the light 
curves are listed in Table 1, which are divided up into seven 


13 The version of Phoebe used for the CV plotting is 2.0a2. 
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Table 1 
All Parameters Used to Calculate Light Curves by Using XRBinary 
Parameters Values Statements 
Fixed 
STAR2TILES 5000 tiles covering star surface 
DISKTILES 20,000 tiles covering disk surface 
BANDPASS (A) 4200 ~ 9000 bandpass of K2 light curve 
Hyow 1.1 power index of Hai 
Lwa z white dwarf luminosity 
Twa e. white dwarf temperature 
Rin? Ue inner radius of disk 
Adjustable 
Components 
orb ak mass ratio 
lord Ue orbital inclination 
Mwa Us white dwarf mass 
Tra e. red dwarf temperature 
Rout Ue outer radius of disk 
Hedge e. Haisk at the outer edge 
€ Ue power index of Tua. 
Lao ets disk luminosity 


Hotspot at the edge of disk 


Tes e. temperature 
Gesmid e centering phase 
Ceswidth sen full width 


Hotspot on the surface of disk 


Gssmin e lower limit of boundary in ¢ direction 


Cssráx U upper limit of boundary in direction 

p m e lower limit of boundary in radius 
direction 

Renaz E upper limit of boundary in radius 
direction 

Tio Ue fractional change in Tai 

Note. 


* Rin can be adjustable for IPs due to the disrupted inner disk. 


fixed parameters and 16 adjustable parameters. Our assump- 
tions for these parameters are as follows. 


1. Although the input and output light curves of XRBinary are 
normalized, the typical luminosity of the white dwarf can 
provide an important reference for the system luminosity of 
the output CV model. Since the mass range for isolated 
white dwarfs is 0.3 Mo to 1.4 Mc; and over 30% of white 
dwarfs are centered on ~0.56 Mẹ (Berg et al. 1992; 
Provencal et al. 1998), the white dwarf masses of all four 
objects were searched in this mass range. By using the 
relation Lwa = 47R24 0744, the white dwarf radius Rwa is 
solely dependent on the white dwarf luminosity Lwa and 
temperature Twa. According to the white dwarf mass-radius 
relationship, as shown in Figure | (e.g., Wood & Horne 
1990; Provencal et al. 1998), Lwa can only be determined 
by Twa for a given Mwa. During the iterations of NMfit, Twa 
calculated from the adjustable parameter Lwa is preset to be 
a fixed parameter. 

2. To improve the reliability of the model, a simple disk 
with minimal free parameters consisting of up to two 
hotspots (one at the vertical side of the edge of the disk 
(es) and the other one on the surface of the disk (ss)) is 
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Rwa (solar radius) 
0.003 0.006 0.009 


0.012 0.015 0.018 


Mwa (solar mass) 


10 20 


30 40 50 
Tyg (x 10° K ) 


Figure 1. Colormap of the luminosity of the primary white dwarf corresponding to Mwa and Twa. The dash line represents the typical mass-radius relationship of the 
white dwarf. The data point marked by RZs refer to the model of RZ Leo derived from the light curve with the phase zero at the secondary dip. The three data points 
marked in red by WDO, WD1 and WD2 refer to the models of WD 1144+011 by using the revised model-0, model-1 and model-2, respectively. The left two data 
points marked by TW and J06 refer to the models of TW Vir and J06324-2536, respectively. 


used to describe the low-inclination CV light curves. This 
CV quiescent disk was assumed to normally extend to the 
surface of the white dwarf (i.e., Rin = Rwa). A boundary 
layer was not considered because it dominates the flux at 
ultraviolet wavelengths rather than at optical wave- 
lengths. A hotspot on the vertical side of the edge of 
the disk can be completely described by three para- 
meters: a uniform temperature Tes, the centering phase 
Gesmia and the full width C.swiatn. Considering that the 
other hotspot on the disk surface can be visible and 
modulate the light curves of low-inclination CVs, this 
hotspot is described using five parameters: Cosmin, Gssmax: 
Rssmin: Rssmax and Tatio. The first four parameters are the 
boundaries of the hotspot, which are the angles and radii 
over which the hotspot extends, respectively. The last 
parameter Tatio is a fractional change in T7; (ie., 
Tatio = ss/Taisk). 

. A different temperature distribution may exist in different 
quiescent accretion disks. The temperature distribution is 
assumed to be a power law in disk radius Raisk; i.e., 
Tüisk X Ri... For a steady-state disk, typically £ = —0.75 
(Wade & Hubeny 1998). However, subsequent papers 
suggested that accretion disks in CVs have much flatter 
temperature profiles with € > —0.75 (e.g., Marsh 1999; 
Orosz & Wade 2003). According to the disk instability 
model (Osaki & Kato 2013a, 2013b, 2014), the quiescent 
DN disk is an optically thin and cool disk with a flat 
radial temperature profile. Observations of various 
quiescent DN (e.g., Z Cha, Wood & Crawford 1986; 
OY Car, Wood et al. 1989; V2051 Oph, Rutkowski 
et al. 2016) confirm that their disk temperature distribu- 
tions are much flatter than the prediction of Ta cx 
Ra. We initially assume € = —0.15 but allow this 
value to be adjusted by NMfit. 

. The adjustable parameter Ly indicating the disk 
luminosity is only calculated from the temperature 


distribution of the disk. This parameter should satisfy a 
default condition that the temperature at the inner edge of 
disk (i.e., Tin) cannot be much larger than 10^ K because 
the average temperature of observed quiescent disks is 
below 104K, which is consistent with the typical 
temperature range of an accretion disk (from several 
10? to 10^ K) predicted by the limit-cycle oscillation 
between hot and cold states (Lasota 2001). 

5. Since both theories and observations suggest that the 
accretion disk around a white dwarf is geometrically thin 
(e.g., Pringle 1991; Frank et al. 1992; Wade & 
Hubeny 1998; Marsh 1999), the height of the disk, 
Hai, is assumed to obey a power low with a power index 
Hyow = 1.1 for all four objects, but Heage is set to be 
adjustable. 

6. Since the components of the model calculated in 
XRBinary (e.g., two component stars, the disk and 
hotspots on the disk) are independent, we can freely 
construct the necessary components to model the 
observed light curves. Three types of models, as listed in 
Table 2, were attempted to verify the necessity of a disk 
around the white dwarf and the hotspots on the disk. 
Model-0 is a detached binary model without a disk. 
Model-1 contains an accretion disk with a hotspot at the 
edge of the disk. Model-2 adds a second hotspot on the 
disk surface. The relative flux contributions in percentage 
from the different model components to the synthetic 
light curves calculated by XRBinary are shown in 
Figure 2. Their zero points are listed in Table 3. 


5. Irradiation Effect Calculated in NMfit 


The irradiation effect is only calculated by XRBinary when 
setting the parameter IRRADIATION to ON. Since a luminous 
secondary may weaken the illumination from the vicinity of the 
white dwarf and the disk, the parameter IRRADIATION is 
preset to OFF if the summation of the luminosity of the white 
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Table 2 
The Models Used in XRBinary 
Component? Model-0 Model-1 Model-2 
White dwarf 1 1 1 
Red dwarf 1 1 1 
Accretion disk 0 1 1 
Hotspot on the disk surface 0 0 1 


Note. 
* The component included in the model or not is indicated by “1” or “0,” 
respectively. 


dwarf and the disk (1.e., LwvatLa) is smaller than the luminosity 
of the secondary L,g. Although XRBinary can calculate the 
heating effect due to the irradiation, the output parameter 7,4 of 
XRBinary is simply derived from L,g by using the formula 
Lia ~ AxRÀ, o TÅ. This means that T,a represents the temper- 
ature of the secondary averaged over its entire surface, which 
can be significantly higher than the true night side temperature 
of the secondary if irradiation in the systems is significant. 
Here, the night side is the surface of the secondary facing away 
from the WD and disk, and the day side is the surface facing 
the WD and disk. 

Davey & Smith (1992) found that the size of the irradiation 
region is usually as large as the entire day side surface of the 
secondary star. The numerical work by Kirbiyik (1982) 
suggested that the difference in the effective temperature 
between the day and night sides of the secondary is very large. 
Considering that the typical difference may reach close to 
10^K (Warner 2003), irradiation would be easily observed if 
this tremendous disparity of the effective temperature exists. 
Therefore, an estimation of the size of the irradiation region on 
the surface of the secondary can help to understand the flux 
contribution from the irradiation region in both CV and pre-CV 
systems. The effective temperature of the secondary on the 
night side can be specified by the semi-empirical CV donor 
sequence (Knigge 2006; Knigge et al. 2011). The derived T,a 
can be regarded as a lower limit to the effective temperature of 
the irradiation region. 

Due to orbital rotation, our viewing angle of the irradiated 
region can change with the orbital phase. Assuming that the 
irradiation region on the leading side of the star'* is 
axisymmetric with the line between the center points of two 
component stars, the flux ratio should reach a maximum 
between phases 0.5 and 0.0; i.e., Fin = fo.s/fo.o, where fo. and 
fo.o are the irradiation flux at phases 0.5 and 0.0, respectively. 
The details of the calculations of F;,, are described in 
Appendix A. By investigating the variations in Fir along with 
a normalized area of the irradiation region (i.e., Si, = Airr/Astar> 
where Air and Astar are the areas of the irradiation region and 
the whole star, respectively), a relation between Fir and Si, 
based on the physical model derived by XRBinary and NMfit 
can be estimated and used to compare with the observed flux 
ratio F,, between phases 0.5 and 0.0. 


6. Comparisons with Gaia Results 
Since the K2 light curves of the four binary systems are from 
the broad bandpass listed in Table 1, their K2 magnitudes 


14 The secondary is assumed to be spherical to simplify the calculations of the 
irradiation effect. 
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(hereafter Kp) translated by Dai et al. (2016) are commonly 
different from their apparent visual magnitudes. Assuming this 
magnitude difference to be a systematic error, the calculated 
magnitudes of the four systems (hereafter CKp2) based on the 
]uminosities of the models can be set to be a reference parameter 
for the comparisons with the results from the Gaia mission (Gaia 
Collaboration et al. 2016). Table 4 lists the distances of all four 
systems based on the absolute stellar parallax in the Gaia database 
(Gaia Collaboration et al. 2018).^ CKp» can be estimated from 
the Gaia distance D, derived from the Gaia parallax (Luri et al. 
2018) and the system luminosity Lan, which is summed for all 
model components and corresponds to the zero point of the 
normalized K2 light curves, by using the following formula, 


D 
CKp, = 251eg| “| BC, 4 siol co (D 
all 


where Lboio = 3.0128 x 10?? erg s L^ and BC, is a bolo- 
metric correction of the main sequence star in the typical 
Johnson-Cousins V band with the solar log g, [Fe/H] and 
[a/Fe]. Note that BC, corresponding to the derived 7,4 can 
cause an uncertainty of CKp for CV and pre-CV systems 
because BC, is improperly used for the white dwarf and disk. 
This calculated uncertainty may explain the discrepancy 
between CKp» and Kp». 

Considering that BC, is a model-dependent quantity with 
many observational constraints, several numerous tabulations 
are provided in the literature. All three popular BC, tables, 
respectively, proposed by Flower (1996), Bessell et al. (1998) 
and Casagrande & VandenBerg (2014) were used to calculate 
CKp» to minimize our errors. The relations of BC, and Tete 
plotted in Figure 3 clearly show that the three BC, tabulations 
significantly differ for cool stars with T.¢ < 4000 K. Inspection 
of Figure 3 indicates that BC, of RZLeo, TW Vir and 
WD 11444-011 are obviously different for each table. Conse- 
quently, the difference between the three BC, tabulations for 
the same Tarp is regarded as the uncertainty of BC,. Compared 
with the two early BC, tables of Flower (1996) and Bessell 
et al. (1998), the updated BC, given by Casagrande & 
VandenBerg (2014) is moderate for cool stars. Therefore, they 
were used to calculate CKp» for all four systems to compare 
with the K2 magnitudes at the zero point of the normalized K2 
light curves (Dai et al. 2016). The interpolated BC, and the 
errors used to calculate CKp; for all four systems are listed in 
Table 4. 


7. Modeling Results and Discussions 
7.1. J06324-2536 


This is a poorly studied DN that experienced several DN 
outbursts in 2009 and 2012 (Korotkiy & Sokolovsky 2012; 
Masi 2012; Ohshima 2012). Dai et al. (2016) showed its double- 
hump light curve extracted from the K2 data archive and several 
quiescent spectra with strong double peaked Balmer emission 
lines obtained from the Large Binocular Telescope (LBT). Based 
on these spectra and the phased K2 light curve, the secondary of 
J0632--2536 is thought to be a K5V star and the orbital 
inclination cannot be lower than 50?. According to the MK 


'S https: / [gea.esac.esa.int/archive/documentation /GDR2 / Miscellaneous / 
sec credit and citation instructions/ 


16 A Zero point luminosity corresponds to an absolute bolometric magnitude 
scale (i.e., Mi.) = 0) recommended by the IAU 2015 Resolution B2. 
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Figure 2. The phased and binned light curves of J06324-2536, RZ Leo and TW Vir superimposed with their best-fitting light curves are plotted in the three left-hand 
panels from top to bottom, respectively. The fluxes of all of the light curves are normalized. The relative flux contributions from different model component for three 
systems are plotted in percentages. The dotted and dashed lines refer to the contributions from two stellar components (white dwarf and red dwarf) and the accretion 
disk, respectively. The solid and dash dot lines denote the contributions from the hotspots at the edge of the disk and on the surface of the disk, respectively. They 
indicate what component is actually contributing to the actual light curve. The three right-hand panels show their corresponding 2D binary configurations at phase 0.75 
by using the Phoebe 2.0. The color denotes the effective temperature. To visualize the hotspot at the edge of the disk, they are filled with black rather than the color 
picked from the color bar because the small temperature difference between the hotspot and the neighboring region of the disk can seriously reduce the contrast of the 


hotspot. The arrow denotes a clockwise rotation direction of the binary system. 


spectral classes (Cox 2000), the initial temperature and mass of 
the secondary are set to be 4410 K and 0.67 Mo, respectively. 
Several high-precision measurements of white dwarf masses in 
CVs (Knigge 2006; Knigge et al. 2011; Zorotovic et al. 2011) 
have indicated that the mean CV white dwarf mass is usually 
larger than that of isolated white dwarfs (Berg et al. 1992; Kepler 
et al. 2007). Thus, for J06324-2536, which lacks an accurate 
measurement of the white dwarf mass, a mean CV white 
dwarf mass of 0.83 Mẹ was preset to be the initial parameter 


(i.e., Mya = 0.83 Mo). Finally, Urban & Sion (2006) used JUE 
spectra of 53 quiescent DN to show that the temperature of a 
white dwarf, Twa, in a CV with an orbital period above the 
period gap is approximately 25,793 K (Sion 1999; Urban & 
Sion 2006). We assume this value for the temperature of the 
white dwarf in J0632--2536 because its orbital period of 
0.314478 day (Dai et al. 2016) is above the gap. 

A large number of trials calculated using XRBinary and 
NMfit suggest that a convergent solution using a normal CV 
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Table 3 
The Zero Points of the Relative Flux Contributions from Different Component 


Components? J0632--2536 RZ Leo TW Vir WD 11444-011 
Binary” 65.1 66.5 20.4 94.3 
Starspot® EE e e 1.1 
Disk" 11.7 0.0 63.0 EE 
Hotspot? 0.0 0.0 0.0 

Hotspot? e 2.6 3.4 

Notes. 


* The relative flux contributions are in percentage. 
b Only consist of white dwarf and red dwarf. 
* Starspots on the red dwarf. 

Accretion disk without hotspot. 


configuration cannot be achieved when setting phase zero to 
the primary dip of the light curve. Hence, phase zero was set 
for the secondary dip. Model-1 was used to model the phased 
light curve of J06324-2536. At first, four parameters (Mya, Twa, 
mass ratio qob» and Tj4 were fixed and the other nine 
parameters consisting of the orbital inclination io» and eight 
parameters of the accretion disk were set to be adjustable. By 
using NMfit, a preliminary accretion disk model indicated that 
the inner radius of the disk is basically close to the white dwarf 
radius, so we then fixed the parameters € and R;,, to derive the 
uncertainties of the other 11 parameters. Since € is an 
insensitive parameter, its uncertainty cannot be obtained. By 
using Phoebe 2.0, a 2D CV configuration at phase 0.75 shown 
in Figure 2 visually indicates that J0632--2536 has a large and 
thick accretion disk with a small hotspot at the edge of the disk. 
Since the orbital inclination of J0632--2536 is not high enough 
for the white dwarf to be occulted by the secondary, the flux 
contribution from the white dwarf is constant and does not vary 
with the orbital period, so this constant flux contribution is 
simply added to the ellipsoidal modulations caused by the 
secondary. Compared with the small white dwarf, the top panel 
of Figure 2 clearly shows that the large and thick accretion disk 
can be partially eclipsed around phase zero. Moreover, the 
hotspot at the edge of the disk, which has a temperature of 6200 
(+100) K and has a phase width of 0.033(2) at phase 0.844(2), 
only contributes a maximum of 5% of the flux from the whole 
system. Since CKp of J0632--2536 listed in Table 4 is almost 
equal to the corresponding Kp», the obtained physical model of 
J06324-2536 is compatible with the Gaia distance. However, 
the Ter of J0632--2536 shown in the Gaia database is 
obviously higher than that derived from the LBT spectra (Dai 
et al. 2016) and the derived 7,4 listed in Table 5. 

The region where the mass transfer stream intersects the 
accretion disk (i.e., the hotspot at the edge of the disk) is simply 
assumed to be a rectangle with a uniform temperature Tes by 
XRBinary. A mass transfer rate Mua (hereafter, the parameters 
with the subscript rd related to the secondary star) can be 
roughly estimated by using the following formula, 


Lace Rout 
G Mwa ' 


Ma ~ (2) 
where Lace is the luminosity of the hotspot at the edge of the 
disk. The mass transfer rate can be estimated to be 
Ma ~ 2.5(+0.2) x 1071? M. yr-!, corresponding to a mass 
loss timescale (i.e., rjj ~ M,a /Ma) of 2.8 x 10? year. The top 
right-hand panel of Figure 4 indicates that M,a = 0.7(1) Mo is 


Dai et al. 


closer to the semi-empirical mass-period relation of Warner 
(2003) than that of Smith & Dhillon (1998). Furthermore, the 
secondary radius R,g = 0.81 is consistent with the radius- 
period relations shown in the bottom right-hand panel of 
Figure 4. The thermal (or Kelvin-Helmholtz) timescale of the 
secondary (i.e., Tkh ~ GMÀ /(LyaRra)) can be calculated to be 
7.6 x 107 year, as listed in Table 6, which is smaller than the 
derived 7. This means that the mass transfer via the L1 point 
is slow and the secondary is always able to maintain thermal 
equilibrium. Therefore, the secondary of J06324-2536 is almost 
indistinguishable from an isolated main sequence star of the 
same mass. The derived mass and radius shown in the bottom 
left-hand panel of Figure 4 confirm that the secondary of J0632 
4-2536 is a typical K5V star (Dai et al. 2016). Inspection of the 
top panel of Figure 2 indicates that the ellipsoidal modulations 
of this K5 star dominate the observed double-hump modula- 
tions. Although the top left-hand panel of Figure 4 indicates 
that 7,4 = 4540(80) K is nearly consistent with a normal K5V 
star based on the MK spectral classes (Cox 2000), this derived 
Tia is around 350 K higher than the prediction of the semi- 
empirical CV donor sequence (Knigge 2006; Knigge et al. 
2011). The most likely explanation for this discrepancy may be 
due to irradiation of the K5V star. Since the irradiation effect is 
not calculated by XRBinary as La > Lwa + La, we investi- 
gated its strength in more detail. Based on our best-fitting 
model, the maximal F;, is around 1.08, as shown in Figure 5. 
As long as 2096 of the irradiation region can be occulted by the 
disk at phase 0.5, the observed flux ratio of J06324-2536 (i.e., 
Fs ~ 0.87) can be explained. Due to our derived large disk 
and moderate orbital inclination, a partial eclipse of the 
irradiation region at phase 0.5 can be expected. 


7.2. RZ Leo 


Dai et al. (2016) confirmed that RZ Leo is a short-period DN 
with an orbital period of 0.07603 day, similar to that derived by 
Patterson et al. (2003) and Kato et al. (2009). Ishioka et al. 
(2001) and Mennickent et al. (1999) indicated that the 
secondary of RZ Leo should be a normal red dwarf, possibly 
a main sequence MO star. Although this spectral type cannot be 
supported by the mass ratio of 0.14 derived from its superhump 
excess of 0.033 (Ishioka et al. 2001), we speculated that 
RZ Leo should consist of a massive white dwarf and a normal 
red dwarf with a small mass of «0.2 Mọ. Thus, we assumed 
that the initial M,4, used in NMfit is 0.15 Mẹ, which means 
Mwa = 1.1 Mo and Tj, ~ 3500 K. Szkody et al. (2017) and 
Dai et al. (2016) detected a short white dwarf spin period of 
~220 s classifying RZ Leo as a member of the IP subtype of 
CVs. From an Ultraviolet (UV) study, Pala et al. (2017) 
determined that the white dwarf had a temperature of 15,014 
(+638) K and contributed 83% of the UV flux. Since the 
humps of the model-1 light curve are located exactly at phases 
0.25/0.75, model-2 with an extra hotspot on the surface of disk 
can be used to describe the offset secondary hump of RZ Leo. 
For the phased light curve with phase zero at the primary 
minimum, a convergent solution was found. However, 
inspection of the bottom left-hand panel of Figure 4 shows 
that the derived secondary mass and radius obviously deviate 
from the predicted mass-radius relation of the semi-empirical 
CV donor sequence (Knigge 2006; Knigge et al. 2011). Based 
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Table 4 
The Results of Four Systems Estimated From the Gaia Database 
Name Test Parallax D, Lan BC, CKp.” Kp 
K mas pc x 10°? erg s^! mag mag 
J06324-2536 5255.5 2.24(5) 446(+10) 13.2 —0.53(4) 14.67(7) 14.4 
RZ Leo Ut 3.6(3) 278(+23) 0.2 —1.4(4) 19.1(4) 17.6 
TW Vir 4850 2.3(1) 435(+19) 5.8 —1.1(1) 16.0(2) 15.9 
WD1144+011 3657.5 4.7(1) 213(+5) 4.8 —1.8(7) 15.5(7) 16.0 
Notes. 


* Interpolated bolometric corrections based on BC, tabulations by Casagrande & VandenBerg (2014). 


Calculated K2 magnitude by using the distance and Ly. 
© Observed K2 magnitude corresponding to Lan- 


on the MK spectral classes (Cox 2000), the derived radius is 
significantly smaller than that of a normal isolated main 
sequence star. Furthermore, the mass and radius also conflict 
with the semi-empirical mass-period and radius-period rela- 
tions of CV secondary (Smith & Dhillon 1998; Warner 2003) 
shown in the two right-hand panels of Figure 4. 

Instead, we created a model with phase zero at the secondary 
dip of the light curve. Initially, we fixed the inner accretion disk 
radius to the radius of the white dwarf. However, a convergent 
solution could not be found. Instead, an accretion disk model 
with two hotspots and Rin much larger than Rwa was found to fit 
the data reasonable well. Ri, > Rwa can be interpreted as a 
truncation of the accretion disk far from the WD, consistent 
with the IP classification of RZ Leo by Szkody et al. (2017) and 
Dai et al. (2016). This means that the physical structure of 
RZLeo may be more complex than the simple disk model 
calculated by XRBinary because accretion curtains near the 
magnetic poles of the white dwarf may be involved. This may 
be the reason that the best-fitting light curve derived by 
XRBinary cannot perfectly describe the secondary hump of 
RZ Leo. Based on this disk model, CKp» calculated by using 
Equation (1) is around 1.5mag larger than Kp). This 
discrepancy is likely to be caused by the lack of an accretion 
curtain in the disk model. In the future, a complete IP model 
including the light from an accretion curtain should help correct 
this difference. In spite of the deviation in the secondary hump 
and CKp;, limits to two key parameters of the disk (Rin and £) 
can be obtained from the double-humped K2 phased light 
curve. 

According to the typical disk-field interaction models (e.g., 
Pringle & Rees 1972; Ghosh & Lamb 1978; White & 
Stella 1988), the corotation radius of a magnetic white dwarf, 
Reo, can be calculated using, 

Reo = f My4 Prot » (3) 
where fo = (GM; ip 41?RÀ)/ 3, R.o and M are in solar units, 
Pot is the spin period of the magnetic white dwarf in unit of 
seconds. By using the parameters listed in Table 5 and a spin 
period of 220s (Dai et al. 2016; Szkody et al. 2017), Reo of 
RZ Leo can be estimated to be 0.085 Ro, which is 17 times 
larger than Rwa = 0.005 Ro. Considering that the accretion 
process in RZ Leo is steady (Dai et al. 2017), Rin is required to 
be almost equal to R.o. However, Rin = 0.211(2) Ro is much 
larger than R.o. In a typical IP (Ghosh & Lamb 1978), Reo is 
regarded as the inner radius of an unperturbed disk and the 
transition region of a truncated disk (i.e., located between R,, 


and Rin) may be actually invisible in the optical band due to a 
possible disruption of the accretion flow caused by the 
magnetosphere of the magnetic white dwarf. Since a hotspot 
on the disk surface always shows R4, < Rin during the 
iterations of NMfit, the parameter of R,smin is fixed to be equal 
to Rin. The high effective temperature of this hotspot indicated 
by Tyatio = 7.2 indicates that more than 97% of the luminosity 
of the disk is from this slim and extended hotspot. Moreover, 
this hotspot is regarded to be a plausible second impact region 
of an inward and overflowing stream striking onto the 
magnetosphere of the magnetic white dwarf. This region was 
previously reported in several IP systems (e.g., EX Hya in 
outburst, Hellier et al. 1989; and QZ Vir in quiescence, Shafter 
& Szkody 1984). Note that all three systems are for an 
unusually short period (under the gap) IPs. Compared to this 
hotspot, the other hotspot at the edge of the disk is much 
smaller and cooler. Both edge and surface hotspots are located 
at the phases 0.723(1) and 0.47(2), respectively. This disk 
configuration derived in quiescence is consistent with the disk 
of EX Hya during outburst, with two hotspots at similar phases 
(Hellier et al. 1989). This disk configuration may be common 
for IP systems. A low-luminosity truncated disk with L; — 
2.3 x 1090 erg s-! (i.e., the luminosity contribution of the disk 
is less than 10% of the whole CV system) may be a 
straightforward conclusion for an IP system like RZ Leo due 
to the lack of a hot inner part of the disk. Compared with the 
other four sources of light in RZ Leo (i.e., the two component 
stars and two hotspots on the disk), the relative flux 
contributions from this truncated disk shown in Figure 2 are 
almost negligible (close to zero). 

The best-fitting CV model indicates that RZ Leo has a moderate 
orbital inclination of 61°0(+0°9) and contains a massive white 
dwarf and a low mass red dwarf. The mass ratio, qo» = 0.080(6) 
is within the error bar of that estimated from its superhump excess 
(Ishioka et al. 2001). The secondary mass and radius (i.e., 
M,a = 0.101(8) Me and Ra = 0.16 Ro) are not only in accord 
with an isolated red dwarf later than M5V but also support the 
semi-empirical mass-period and radius-period relations of CV 
secondaries (Smith & Dhillon 1998; Warner 2003). The small 
mass transfer rate, Ma ^ 3.01 (+0.07) x 107? Mo yr-! calcu- 
lated by using Equation (2), seems to explain the small and cool 
hotspot at the edge of the disk that is shown in the disk model. 
Recently, Dubus et al. (2018) showed the average mass transfer 
rates of a sample of ~130 CVs with a parallax distance in the 
Gaia DR2 catalog. The estimated mass transfer rate of RZ Leo lies 
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Figure 3. Diagram of bolometric correction in V band (BC,) against the effective temperature of a main sequence star (Tete). Three legends of Casagrande14, 
Bessell98 and Flower96 refer to the BC, tabulations given by Casagrande & VandenBerg (2014), Bessell et al. (1998) and Flower (1996), respectively. The four 
dotted lines from left to right denote 7;4 of WD 1144+011, RZ Leo, TW Vir and J06324-2536, respectively. 


well within the region where a stable, cool disk exists within the 
system. This explains a lack of DN outbursts for RZ Leo. 
Tar ~ 3.4 x 10!° year calculated from the mass transfer rate is 
larger than Tkn of the secondary. The derived secondary mass and 
radius further confirm that the secondary of RZ Leo is in thermal 
equilibrium, since they are consistent with the mass-radius 
relation of the isolated main sequence stars shown in the bottom 
left-hand panel of Figure 4. Since 7,4 = 3710(430) K suggests a 
slightly earlier spectral type than MSV, this higher Tj, may be 
attributed to irradiation of the companion. Due to Lig > Lwa + La, 
the irradiation effect is not calculated by XRBinary. However, the 
maximal Fj, ~ 1.54 that is shown in Figure 5 indicates a 
significant irradiation effect in RZ Leo. A prominent irradiation 
effect can be expected in an IP system because the white dwarf 
primary is not embedded in the truncated thin disk. Furthermore, 
the small disk of RZ Leo implies that the irradiation region is 
visible at phase 0.5 despite the similar orbital inclination to that of 
J06324-2536. Although Fəbs ~ 0.93 does not seem to support a 
large irradiation effect, the presence of an accretion curtain and a 
brighter second impact region complicate an accurate calculation 
of the irradiation effect in RZ Leo. 


7.3. TW Vir 


By using infrared photometry, the secondary is estimated to be 
a M3V star with a mass of 0.43 Mọ and a radius of 0.48 Ro 
(Mateo et al. 1985). According to the MK spectral classes 
(Cox 2000), the initial temperature of the secondary was set to be 
3367 K. Based on the mass ratio of 0.44 (Shafter 1983), the white 
dwarf mass of TW Vir is calculated to be 0.98 Mo. Since no 
accurate white dwarf temperature has been obtained from the UV 
spectra (Córdova & Mason 1982; Szkody 1985; Hamilton 
et al. 2007), and TW Vir is a long period (above the gap) system 
with an orbital period of 0.182682(3) day derived from the K2 
data Dai et al. (2017), we set the same initial white dwarf 


temperature as that of J0632--2536. Because the secondary hump 
of TW Vir is not located at orbital phase 0.25 or 0.75, model-2 is 
used for the calculations. Two phased light curves corresponding 
to phase zero at the primary and at the secondary dips were used 
to derive two disk models, respectively. In these two disk models, 
three parameters involving the positions of two hotspots on the 
disk in Ç direction (i.e., Gesnia, Gssmin aNd Cosmax) Show nearly 
constant phase differences of 0.5, while all of the remaining 
parameters are almost the same. Therefore, both models are not 
significantly different and the double-hump modulation of 
TW Vir is mainly caused by the two hotspots on the disk. 
Since the xX? of the model based on the phased light curve 
with phase zero around the secondary dip is slightly smaller 
than that with phase zero at the primary dip,the primary hump 
was set to phase at 0.75 and the corresponding model is 
analyzed in the following discussion. The two best-fitting 
orbital parameters:qoy = 0.41(40.05) and i, = 44°3(+0°5) 
listed in Table 5 are consistent with the previous results derived 
from the infrared and UV data (Córdova & Mason 1982; Mateo 
et al. 1985). The derived mass ratio confirms the previous 
classification of TW Vir as a U Gem type DN (O’Connell 1932) 
but the derived orbital inclination is smaller than that obtained 
by Hamilton et al. (2007). The calculated CKp» listed in 
Table 4 suggests that the best-fitting model of TW Vir is 
consistent with the Gaia distance. The relative flux contribu- 
tions from the hotspot at the edge of the disk around phase 0.75 
shown in Figure 2 almost perfectly reproduces the significant 
primary hump, which lasts for over half of an orbit. Although 
the other hotspot appearing on the outer radius of the disk 
covers a long phase width of 0.42,the small Tio = 1.20 
(+£0.02) may explain its small relative flux contributions, as 
shown in Figure 2. Like RZ Leo,the relative flux contribution 
from the disk is close to zero. The derived mass and radius of 
the secondary fits with the semi-empirical secondary mass- 
period and radius-period relations (Smith & Dhillon 1998; 
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Table 5 

Photometric Solutions for Three DN 
Parameter* J0632+2536 RZ Leo TW Vir 
Model type” Model-1 Model-2 Model-2'* 
Orbit 
Qon (Ma/ Mya) 0.8(1) 0.080(6) 0.41(5) 
iorb (degree) 60.3(5) 61.0(9) 44.3(5) 
White dwarf 
Mya(M 5) 0.81(8) 1.26(3) 1.10(3) 
Rya( R5) 0.01 0.005 0.007 
TG 254 15.0 19.1 
Lua 15 0.045 0.23 
Red dwarf 
MaMa) 0.7(1) 0.101(8) 0.45(6) 
Rra(Ro)° 0.81 0.16 0.48 
Tra 4.54(8) 3.71(3) 4.00(4) 
Lu 9.6 0.17 2.0 
Accretion disk 
Rin(Ro) 0.01? 0.211) 0.007* 
Rou(Ro) 1.22(1) 0.401(2) 0.536(3) 
Hedge(Ro) 0.091(4) 0.0801(6) 0.1174(3) 
E —0.11* —0.12(6) —0.11(4) 
Lao 3.9(3) 0.0006 1(6) 3.36(6) 
La 4.1 0.023 4.0 
Hotspot at the edge of the disk 
Tas 6.2(1) 3.71(2) 4.80(2) 
Cesmia (phase) 0.844(2) 0.723(1) 0.727(6) 
Ceswidth (phase) 0.033(2) 0.057(4) 0.296(3) 
Hotspot on the surface of disk 
Coginiti (phase) 0.38(2) 0.32(5) 
Cssmax (phase) 0.56(2) 0.74(3) 
Rssmin(Ro) 0.211° 0.489(8) 
Rysmax(Ro) 0.217(5) 0.536° 
Tio 7.22) 1.202) 
x? 11.9 2.3 0.66 
Notes. 


? The unit of temperature and luminosity is 10° K and 10? erg s^ 


1 


, respectively. 


5 The model with a superscript “irr” denotes that the irradiation effect is included in XRBinary. 


* Calculated by XRBinary. 
4 Fixed in NMfit program. 


© Insensitive to the observed light curves. 


Warner 2003). This indicates that the spectral type of the 
secondary is close to M2V, which is similar to the result from 
infrared photometry (Mateo et al. 1985). However,compared 
with the K5V star of J06324-2536 and the M5V star of RZ Leo, 
the M2V star of TW Vir never shows large-amplitude 
ellipsoidal modulations. Instead, it gives rise to a small 
amplitude second hump in the phased light curve. By using 
Equation (2), the estimated Ma listed in Table 6 indicates that 
Tm ^ 9Ì Tkn, Which is common for most CV secondaries 
(Patterson 1984). A somewhat oversized secondary of TW Vir 
relative to an isolated main sequence star of the same mass is 
not obtained, implying that M,a may be overestimated. 

Like J06324-2536, T.i = 4850 K listed in Table 4 is higher 
than the derived 7,4, = 4000(+40) K. This may imply that the 


10 


higher Tore derived by Gaia is common for DN due to the 
possible contributions from hotter components in these systems 
(e.g., the white dwarf and accretion disk). Since 7;4 is obviously 
higher than the effective temperature of the M2V star shown in 
Figure 4, we investigated the irradiation effect of the secondary 
in this system. Note that TW Vir is the only object with 
Lia < Lwa + Lg. Accordingly, the calculation of irradiation was 
included in XRBinary. However, both light curves calculated 
with and without irradiation by XRBinary are almost identical. 
Hence, irradiation in TW Vir may be weak, which is also 
demonstrated by Fp; ~ 1 with two minima in the light curve at 
almost the same flux level. Since Lq calculated from a near-flat 
temperature distribution with a power index of —0.11(-£0.04) is 
nearly two times L,g, and the total contributions of two hotspots 


THE ASTRONOMICAL JOURNAL, 156:153 (17pp), 2018 October Dai et al. 
_ 
zi : - 
7.0 -~ p 0.0 

pon wb2 m 

— 9s qe Wb1 & 

x 6.0 AT o 
e i Ty oo TOS z 

o 4 IU a 
ro) —wp2— 7 m 
= 5.0 H eL F 
x " yer o 
— P 7 -0.6 = 

P PP 
= 40 E 3 
fd o 
3.0 - i pm 
"U[. p --- MK spectral classes |; RÉS m Warner (1995) 
CV donor sequence |, — - Smith & Dhillon (1998) 
29S 
"WB 4 0.0 

A^ T 
E 5 

E: S 
= -0.3 5 

E Qa 
[e] un 
2 $ 
p S 
zA 06 g 
= 2 
D = 
2 2 


-0.9 


-0.6 -0.3 
Logio My (solar mass) 


Warner (1995) 
— - Smith & Dhillon (1998) 4.0.9 


2 4 6 8 
Orbital Period (hr) 


10 


Figure 4. Four relationships of the secondaries. Top left-hand panel: the relationship of the logarithm of mass and the effective temperature. Bottom left-hand 
panel: the logarithm of mass-radius relationship. Top right-hand panel: the period—mass relationship. Bottom right-hand panel: the period-radius relationship. The 
dashed and solid lines denote the relationships based on the isolated low-mass stars (Cox 2000), and the semi-empirical CV donor sequence (Knigge 2006; Knigge 
et al. 2011), respectively. The dashed—dotted and dotted lines describe the relationships derived by Smith & Dhillon (1998) and Warner (2003), respectively. The 


symbols of the data points are the same as those used in Figure 1. 


Table 6 
The Estimated Parameters of the Secondaries of Four Systems 
Name? Ma TM Tkn TM / Tkh 
10719 Mo yr! 10? year 10? year 
J06324-2536 2.5(£0.2) 2.8(+0.5) 0.076 37 
RZ Leo 0.0311 32.4(£2.7) 0.45 72 
(+0.0008) 

TW Vir 1.92(0.05) 2.3(+£0.3) 0.25 9 
WD 1144 67.8(::16.0) 0.1(2-0.03) 0.006 18 

+011° 
WD 1144 7.9(£1.8) 1.0(+0.3) 0.027 37 

+011 
Notes. 


* Calculations by using the parameters listed in Tables 5 and 7. 
b 

For the model-1. 
* For the model-2. 


are only 16% of Ly, the weak irradiation effect of the secondary 
is overpowered in the light curve. Large variations in the disk 
luminosity may further weaken the ability to detect irradiation. 
There are, however, two issues which further complicate our 
conclusions. The first is that the orbital modulation of TW Vir 
shown in Figure 2 is only extracted from the part of the quiescent 
light curve around its superoutburst. The second is that the 
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double-hump modulation does not always maintain stability, and 
the secondary hump has a small amplitude and shows a large 
scatter in the unbinned light curve (Dai et al. 2017). Further 
analysis of its light curves at different times may reveal more 
details of the irradiation effect in TW Vir. 


7.4. WD 11441011 


This poorly understood variable star is classified as a 
DA+dMe binary from a single optical spectrum obtained by 
Berg et al. (1992). For comparison with this spectrum obtained 
over 20 years ago, a new spectrum was taken on 2017 January 
22 by using the BFOSC spectrograph attached to the Xinglong 
Observatory 2.16 m telescope (XL 216, Fan et al. 2016). This is 
a better optical spectrum that shows some changes from the 
original. The details concerning our spectrum are presented in 
Appendix B. A synthetic analysis based on the single-hump 
light curve extracted by Dai et al. (2016) can test the possibility 
of the existence of a disk around the primary white dwarf of 
WD 1144+011. Considering that WD 1144+011 is a long 
orbital period system with an orbital period of 9.81 hr (Dai 
et al. 2016), the initial parameters of the white dwarf are set to 
be the same as those of J0632+2536. Combined with the 
parameters of a M dwarf, all three models were used to attempt 
to reproduce the stable single-hump modulation that is apparent 
in the K2 light curve. 
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Figure 5. Diagram of the flux ratio between phases 0.5 and 0.0, Fir, changing along with a normalized area of the irradiation region Sir. The red symbol WD denotes 
Sir ~ 0.05 corresponding to F,,, = 1.03 estimated from the single-hump modulation of WD 1144+011. 


7.4.1. Model Without a Disk 


In principle, model-0 can only produce ellipsoidal modulation 
caused by the orbital motion of the Roche-lobe filling secondary. 
The asymmetry shown in the single-hump modulation of 
WD 1144+011 cannot be explained by a pure ellipsoidal 
modulation derived from model-0. Moreover, the irradiation 
effect in WD 1144+011 that was calculated by XRBinary cannot 
resolve this problem because its long orbital period may imply a 
bright secondary like J0632--2536 (i.e., Liq > Lwa). In spite of 
this, an irradiation region on the surface of the secondary was 
added to allow for an asymmetrical single-hump modulation (i.e., 
a bump on the rise to the maximum). We assumed that this 
irradiation region is a circular bright starspot on the secondary, 
which can be easily calculated by XRBinary. Since the single- 
hump modulation of WD 1144+011 observed in K2-C1 is very 
stable (lasting at least a complete campaign periods of ~3 months, 
Dai et al. 2016), a stable region due to irradiation of the secondary 
may be more plausible to explain this long-term steady 
modulation than a variable starspot on an active M type dwarf. 

In XRBinary, a starspot on the secondary can be described by 
the four parameters (i.e., 0P, ¢*?, RP and T®.,)'’. For WD 1144 
+011, this starspot is assumed to be a hotspot (i.e., T, P. 1). 
All of the parameters including four component parameters 
(Mwa, Twa» orb and 7,4) and four parameters of the starspot were 
set to be adjustable in NMfit. The derived Mwa = 1.11(8)Mz is 
almost the same as that of TW Vir shown in Figure 1, while 
Twa = 27,100 K suggests a hotter and brighter white dwarf. 
Compared with the three DN, WD 1144+011 is the only object 
with a slightly cooler M dwarf (T;4 = 3500(+100) K) than the 
predicted mass-temperature relation of the semi-empirical CV 
donor sequence (Knigge 2006; Knigge et al. 2011). This is 
consistent with Tog = 3657.5 K shown in the Gaia database. 
Compared with the higher Tz of the two DN J06324-2536 and 
TW Vir, this compatible 7.; can be explained by the lack of a 
disk around the white dwarf. Since the companion star in 


17 OP ang oP are the coordinates of the spot center in the spherical polar 
coordinate system. ¢°? = 0 direction lies in the orbital plane and points in the 
direction of motion of the secondary in its orbit. à? increases in a right-handed 
sense about the 0"? = 0 direction. R® is the angular radius of the spot as seen 
from the center of the secondary star. T% is the ratio of the starspot 
temperature to the local effective temperature of the unspotted star. 
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WD 1144+011 is cooler than the secondaries in the other 
systems discussed here, it has the largest uncertainty in BC,, as 
shown in Figure 3. This is reflected in the large range of derived 
CKp; (14.8 ~ 16.2 mag), which is consistent with the corresp- 
onding Kp; value listed in Table 4. Like TW Vir, the ellipsoidal 
modulations caused by the M star never dominates the orbital 
modulation. The mass and radius of the secondary deviate from 
all the mass-radius, mass-period and radius-period relations 
shown in Figure 4. Note that both semi-empirical mass-period 
and radius-period relations are only available in the period range 
of 1.3 ~ 9 hr (Smith & Dhillon 1998; Warner 2003). With an 
orbital period of 9.81 hr, WD 1144+011 is beyond this period 
range. The bottom left-hand panel of Figure 4 shows that the 
inflation of the secondary is about 38%. Besides the known 
inefficient mechanisms of donor bloating (e.g., tidal deformation, 
rotational deformation and irradiation of the secondary) 
discussed by Knigge et al. (2011), the deviation from thermal 
equilibrium of the secondary caused by the donor mass loss can 
result in a large donor inflation of up to 20% ~ 30% (Patterson 
2005; Knigge 2006). However, this mass loss mechanism is still 
not enough to explain the derived inflated secondary star. 
Although both light curves calculated with and without 
irradiation by XRBinary are almost identical due to Lia > Lwa, 
the final x? of model-0™ is slightly smaller than that of model- 
0. Thus, the best-fitting parameters and their uncertainties listed 
in Table 7 are calculated using model-0"". The reproduced light 
curves and the corresponding 2D pre-CV configuration are 
shown in the top left-hand and right-hand panels of Figure 6, 
respectively. This detached binary configuration with a low 
inclination of 14°3(+0°4) supports the previous classification 
of WD 1144+011 as a pre-CV system (ie. a DA+dMe 
detached system). Our model also shows that the orbital 
variation seen in the light curve is dominated by the hotspot on 
the secondary's surface. Compared with the three DN, the pre- 
CV WD 1144+011 is the only object with Fobs = 1.03 larger 
than 1. Thus, T;T, = 1.097(£0.003) combined with Fobs 
suggests that the calculated irradiation region may be a good 
representation of the physical picture of WD 1144+011. The 
effective temperature of the starspot can be estimated to be 
around 3900 K based on T; 3... A relation between Fir and Sirr 
shown in Figure 5 indicates weak irradiation in WD 11444-011, 
similar to J06324-2536, and Sir ~ 0.5 corresponding to 
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Figure 6. The phased and binned light curve of WD 1144+011 superimposed with the best-fitting light curves denoting with the red, black and blue lines calculated 
by using the model-0, model-1 and model-2 are plotted in the top, middle and bottom left-hand panels, respectively. The solid and dash lines respectively refer to two 
types of calculations including and excluding the irradiation effect in XRBinary. The three right-hand panels from the top to bottom show the corresponding 2D binary 
configurations at phases 0.6, 0.6 and 0.75, respectively. All symbols are the same as that of Figure 2. 


Fos = 1.03 is consistent with the size of the starspot calculated 
from R® = 27°7(40°7). Due to the lack of an accretion disk, 
this irradiation region is much smaller than the typical size of 
an irradiation region in a CV (Davey & Smith 1992). 


7.4.2. Model With a Disk 


Since an extra light source is imperative to explain the 
asymmetrical single-hump modulation, model-1 and model-2 were 
investigated. All of the adjustable parameters and their uncertain- 
ties are listed in Table 7. Mya is estimated for both models and is 
given in Table 6. The white dwarf mass of model-1 is smaller than 
the prediction of the average white dwarf mass in a CV with the 
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same orbital period (Ritter & Burkert 1986; Zorotovic et al. 2011). 
Twa from both of the disk models is consistent with Twa from the 
non-disk model, suggesting that WD 11444-011 does contain a hot 
white dwarf. The final x of the two disk models are almost the 
same as that of the non-disk model-0. All three models imply an 
oversized secondary with generally consistent masses and radii. 
The top left-hand panel of Figure 4 shows that the secondary of 
model-1 has an extremely high effective temperature, which is 
over twice the average effective temperature of an M dwarf. 
Although the secondary of model-2 agrees with a G9V star based 
on the MK spectral classes (Cox 2000), a spectral type of GOV 
contradicts with our spectrum. Compared with the Gaia distance 
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Table 7 
Three Photometric Solutions for the Pre-CV WD 1144+011 


Parameters? WD 11444-011 
Model type? Model-0"" Model-1 Model-2 
Orbit 
dorb(Mra/Mwa) 0.49(2) 0.99(2) 0.92(9) 
iorb (degree) 14.3(4) 17.7(4) 26.3(6) 
White dwarf 
MyaMo) 1.11(8) 0.68(16) 0.9(2) 
Rya( R3) 0.007 0.012 0.009 
Tia 27.1 27.0 25.7 
wd 0.91 2.5 1.3 
Red dwarf 
Mra(Mo) 0.54(4) 0.67(16) 0.8(2) 
Rya(Ro) 0.87 0.97 1.03 
Tra 3.5(1) 7.44(7) 5.23(6) 
Lia 4.1 100.7 27.7 
Bright starspot on the red dwarf 
0? (degree) 38(3) 
o? (degree) 0.2* 
R* (degree) 27.7(7) 
TP. 1.097(3) 
Accretion disk 
Rin(Ro) 0.01° 0.0094 
Rou(Ro) 0.597(5) 0.808(8) 
Heage(Ro) et 0.363(8) 0.176(3) 
€ e —0.11* —0.18(4) 
Lao e 13.7(42.7) 13.5(42.3) 
d e 23.0 19.4 
Hotspot at the edge of the disk 
Tes 6.50(5) 5.38(7) 
Gesmia (phase) 0.623(4) 0.548(8) 
Ceswidth (phase) 0.412) 0.17(1) 
Hotspot on the surface of disk 
Gesinin (phase) 0.129(6) 
Gssmax (phase) 0.221(2) 
Rssmin(Ro) 0.605(1) 
Rssmax(Ro) 0.808° 
Tatio 1.90(4) 
x 8.6 7.9 8.8 


Note. All footnotes are the same as those of Table 5. 


of WD 1144+011 listed in Table 4, both L,a for the model-1 and 
model-2 are too large. Hence, the consistency of the effective 
temperature may be only coincidental. To recheck the deviations 
in 744, in model-1 and model-2, 7,4 was fixed to be 3500 K in 
line with the secondary temperature found using model-0. 
However, no convergent solution could be obtained. In spite of 
Lia > Lwa + Lg, irradiation was still included in these models. 
Like model-0, the light curves calculated with and without 
irradiation for the two disk models only show small discrepancies 
around phases 0.5 and 0.0. Both calculated F; are far larger than 
Fas. This means that the assumed large irradiation of the 
secondary is not supported by the K2 data. Compared with the 
pre-CV configuration obtained from model-0, the two models with 
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an accretion disk are not convincing. As such, our modeling based 
on the K2 phased light curve provides further evidence for the 
classification of WD 1144--011 as a detached pre-CV system. 


8. Conclusions 
8.1. Synthetic Codes 


Based on the light curve synthesis code XRBinary derived 
by E. L. Robinson, NMfit was developed to analyze the light 
curves of the four low-inclination systems: J0632+2536, 
RZ Leo, TW Vir and WD 1144+011. All of the parameters 
of the best-fitting models and their uncertainties are estimated 
by NMfit. Additionally, Phoebe 2.0 was used to visualize the 
configuration of each system using each systems’ best-fitting 
parameters. Since phase zero is hard to identify in any low- 
inclination CV system with a double-hump modulation due to 
the lack of a significant eclipse feature, we tested models that 
had phase zero located at either the primary or secondary 
minimum, and we chose the model which matched the 
observed light curve best. Except for WD 1144+011, which 
had a single maximum in its phased orbital light curve, the 
derived CV models of the other three DN indicate that phase 
zero should be placed at the secondary minimum. 


8.2. Physical Models 


For J0632--2536 and TW Vir, the best-fit disk models show 
that the primary hump is mostly due to the hotspot at the edge 
of the disk, which is a key indicator of mass transfer via the L1 
point. Another hotspot on the disk surface can explain the 
phase difference between the two humps of the double-hump 
modulations. For WD 1144+011, the bright starspot represent- 
ing irradiation of the secondary star is responsible for the 
modulation seen in the optical light curve. 

The lack of an accretion curtain in the disk model of RZ Leo 
implies that this model may not be appropriate for a 
comparison with the measured distance by Gaia. The derived 
physical models of the other three binary systems are consistent 
with the results from Gaia DR2. The flat power law index of 
the disk found in all three DN models (£ ~ —0.11) is similar to 
previous observations (e.g., Wood & Crawford 1986; Wood 
et al. 1989; Rutkowski et al. 2016) and supports the theory that 
a quiescent CV disk deviates from the temperature distribution 
of a typical steady-state disk. A low-luminosity accretion disk 
model of RZLeo derived from its K2 light curve further 
confirms that RZ Leo is an IP system with two hotspots on a 
truncated disk. One of the hotspots on the disk surface 
contributes a significant fraction of the disk luminosity (79796 
of La), and is located at the inner edge rather than the outer 
edge of the disk. This may be evidence of an impact region 
between an inward and overflowing stream and the magneto- 
sphere of the magnetic white dwarf. Compared with RZ Leo, 
the small hotspot of J0632--2536 and two hotspots of TW Vir 
covering large phase ranges are only small contributors to the 
disk luminosity. Our spectrum of WD 1144+011 with a 
relatively flat continuum and H,, emission supports its previous 
classification as a DA+dMe system (Berg et al. 1992). We 
note that WD 11444-011 shows different flux levels in the 
continuum and emission lines. The model light curve based on 
the asymmetrical single-hump modulation requires an extra 
light source (i.e., a weak irradiation region of the secondary 
rather than a large hotspot at the edge of the disk) to explain the 
modulation of WD 11444-011. 
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6.3. The Secondaries 


The estimated M,4 for all four objects are within a range 
of 107? ~ 107? Mo yr-!. Except for the pre-CV WD 1144 
+011, which contains an oversized secondary, the other three 
DN have secondaries in thermal equilibrium with masses 
and radii conforming to the semi-empirical CV donor sequence 
(Knigge 2006; Knigge et al. 2011) and MK spectral classes (Cox 
2000). The derived effective temperatures of all three DN 
are significantly higher than predicted. Hence, a DN system 
containing a substantially hotter secondary may be a common 
feature rather than a peculiarity. This can be attributed to 
irradiation of the secondary because 7,4 calculated from L4 is an 
average parameter that can be increased by irradiation. Compared 
with Terr listed in the Gaia catalog, the lower 7,4 of the two DN 
J0632--2536 and TW Vir may be due to contamination from a 
hot white dwarf and disk. This is further supported by the 7;4 of 
the pre-CV WD1144--011, which is almost consistent with the 
Gaia Tz. It should also be noted that the Gaia temperatures are 
determined from three broad bandpasses (Andrae et al. 2018) and 
the DR2 release notes urge caution in using them.'? 

Although the double-hump modulation of J0632--2536 can 
be explained by the partial occultation of the irradiation region 
on the surface of the secondary due to a large disk and a 
moderate orbital inclination, investigation of irradiation in the 
other two DN implies that the effect of irradiation in a CV 
system is complicated and blended with other modulations. The 
flux contribution from the secondary of TW Vir is the lowest 
(i.e., Lia < Lwa + La) among all four objects. Weak irradiation 
may exist in the DN TW Vir and the pre-CV WD 11444-011. 
The former can be further tested by additional light curves 
obtained when the double-hump variation is evident, while the 
latter can be further checked by taking a time series of spectra 
over the course of the complete orbital period of 9.81 hr. 
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Appendix A 
A Simple Calculation for the Irradiation Effect of the 
Secondary in a Low-inclination CV System 


Assuming that the secondary is a sphere and the irradiation 
region with a higher effective temperature of Ti, is axisym- 
metric with the line between the center points of the two 
component stars, the area of the complete irradiation region Sirr 
can be calculated by the following formula, 


Ain = 2nRÀ(1 — cos), (4) 


35 https: / [gea.esac.esa.int/archive /documentation/ GDR2/ 
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Figure 7. This sketch displays an irradiation region of the secondary star with a 
segment shape. The center of the secondary is at O point. The OM line is 
located on the orbital plane, while the OS line indicates a plane vertical to the 
line of sight. 0 is the half of the opening angle of the irradiation region. 


where R,q and 0 are the radius of the secondary and the half 
opening angle of the irradiation region shown in Figure 7, 
respectively. Therefore, 0 should be in a range of 0 ~ 7/2. Due 
to the projection effect (i.e., the orbital inclination less than 
90°), the orbital modulation of the secondary will cause the 
fraction of the irradiation region seen by an observer to vary 
over the orbital period. The total flux at an orbital phase @ can 
be described as 


4 
o Tas. 


f =Ain(d)o T + TRÀ, — ALO) 


(5) 


where Aj (@) and Ttar are the visible area of the irradiation 
region at the phase @ and the effective temperature of the 
unirradiated part of the secondary, respectively. Note that 
Equation (5) never considers the occultation of irradiation 
region by the white dwarf or the disk. In principle, A; (0.0) is 
the minimum, while A, (0.5) is the maximum. An inspection 
of Figure 7 indicates that the former can be estimated by the 
formula, 


0 
Ai, (0.0) = aR [f 0.sn)db 02i 
0 0<i 


(6) 


where i and wy are the orbital inclination and the second 
coordinate of the spherical coordinate system, respectively. 
Moreover, 6” satisfies the following relation, 


Rrasiny cos0^ — Ryacosv 
FREI 
Ry cosi 


(7) 


Ry sini 


By combining Equations (6) and (7), AL. (0.0) can be expressed as, 


02i 
0ci 


2R 2 Cu 


Ai, (0.0) = f (8) 
0 
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Figure 8. The optical spectrum of WD 1144+011 covering a wavelength range from 3900 A to 8300 A was taken by using the Xinglong 2.16 m telescope on 2017 


January 22. The characteristic spectral lines are marked. 
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where Cir = I arccos( 7) - siny% dy. Based on A, (0.0), 


Aj.,(0.5) can be expressed as, 


Ain m 
Air 


A1, (0.0) 02i 


0ci 9» 


Ain (0.5) = f 
By using Equation (5), the flux ratio between the phases 0.5 and 
0.0 (Le., Fir = fo.s/fo.o) can be calculated for a given 0. Based on 


the normalized area of the irradiation region, Sir = Ain/Astar (i.e., 
cos 0 = 1 — 25;,,), Fi, can be described as, 


Fi 
(T 27 Sirr CT sar Hae, Ci) yn 0 > i 
(T — Cin) T, du + Cir Tas 
e 2549] =| 20 + Sin) 0«i 
star 
(10) 


The relations of F;,, versus Sir for the three DN and one pre-CV 
system discussed here are plotted in Figure 5. 


Appendix B 
A New Observed Optical Spectrum of WD 1144+011 


A 1.8 arcsec slit and the G4 grating with a low resolution of 
1953 (i.e., 2.97 A pixel ') was used with an exposure time of 
1800s. Flux standards and FeAr lamps were used along with 
IRAF'? reductions to produce the final calibrated spectrum 
shown in Figure 8. 

This spectrum is dominated by a strong and fairly broad Ha 
emission line and several molecular absorption bands (e.g., 
MgH5211 À and CaH6358 A), and is similar to the previous 


19 TRAF is distributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science Foundation. 
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spectrum obtained by Berg et al. (1992). Thus, our spectrum 
supports the previous spectral classification of WD 1144+011 
derived by (Berg et al. 1992). An active M star can show 
variable H,, in emission. However, our spectrum shows a flatter 
blue continuum with a higher red flux level. The measured 
average equivalent width of H, is close to 100 A, although this 
line is difficult to measure accurately due to the uncertainty of 
the continuum. It is clear from both our spectrum and that of 
Berg et al. (1992) that there is no prominent blue continuum 
from a hot white dwarf or disk and that the Balmer lines are in 
emission and of variable strength. The changes in flux levels 
between our spectrum and the spectrum presented in Berg et al. 
(1992) could be due to phase differences when the two spectra 
were taken. 
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